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Abstract 

New approach to highly accurate calculations of the photodetachment cross-sections of the neg¬ 
atively charged hydrogen ions is developed. This effective, fast and numerically stable method is 
based on the use of the Rayleigh’s formula for the spherical Bessel functions. Photodetachment 
cross-sections of the negatively charged hydrogen ion(s) °°H~, 1 H _ (protium), 2 H~ (deuterium D~) 
and 3 H~ (tritium T~) are determined with the use of highly accurate, truly correlated, variational 
wave functions constructed for these ions. Our method allows one to investigate the problem of 
photodetachment of the negatively charged hydrogen ions completely, i.e. without any limitation 
by numerical values of the photo-electron momentum p e . 
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I. INTRODUCTION 


The negatively charged hydrogen ions °°H~ (nuclear mass is infinite), 1 H _ protium) and 
2 H~ (deuterium, or D “) are of great interest in stellar astrophysics 1] - 4]. These ions 
determine the absorption of infrared and visible radiation in photospheres of all stars, if 
temperatures in their photospheres are restricted between T max ~ 8,250 K (late A-stars) 
and T min & 2,750 K (early M-stars). In the late F, G and early K stars the absorption 
of infrared and visible radiation by the negatively charged hydrogen ions is maximal. This 
includes our Sun which is a star of spectral type G2. The great role of the H~ ions for 


0, II 


and 


our Sun was suggested by R. Wild in 1939 (see discussions and references in 
4]). An effective absorbtion of large amount of infrared solar radiation by the negatively 
charged hydrogen ions is crucial to reach a correct thermal balance at our planet. The 
leading contribution to the light absorption by the negatively charged hydrogen ion H _ is 
related to the photodetachment of this ion 

H _ + hv = H(ns) + e~ (1) 

where the notation hv designates the incident light quantum, while e~ means the free elec¬ 
tron, or photo-electron, for short. The notation H(ns) means the hydrogen atom in the 
bound ns—state, where n is the principal quantum number and s means s —states with 
l = 0 where i is the angular momentum of the hydrogen atom. The main goal of this 
study is to determine and investigate the photodetachment cross-sections of the negatively 
charged hydrogen ion(s). A large amount of preliminary work for this problem, including 
derivation of different formulas needed for numerical calculations of the photodetachment 


cross-sections of the H ions can be found in 


[sj . The paper 


also contains a review of the 


photodetachment of negatively charged hydrogen ions and extensive bibliography on this 
subject. This allows us to restrict our analysis in this study to the discussion of computa¬ 
tional method and results obtained with the use of this method. Note also that a very good 
discussion of the photodetachment of the negatively charged hydrogen ions can be found in 
some textbooks, e.g., in 0, & and Q|. 

Our approach developed in this study is based on the use of highly accurate, three-body 
(or two-electron) wave functions constructed for the negatively charged hydrogen ions. Such 
wave functions are truly correlated. The final state wave functions are represented in the 
form of products of two single-electron wave functions, i.e. they are not correlated wave 
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functions. In general, the overlap integrals between truly correlated and non-correlated 
wave functions cannot be computed without a substantial loss of accuracy. For the pho¬ 
todetachment of the negatively charged hydrogen ions this problem, i.e. loss of overall 
accuracy, is known since the first works by Chandrasekhar [q| - [11], In this study we avoid 
any loss of high accuracy during calculations of the three-particle integrals with spherical 
Bessel/Neumann functions by using our recently developed, original method for computa¬ 
tion of such integrals written in relative coordinates. This method allows one to perform 
accurate numerical calculations of the photodetachment cross-sections for different hydrogen 
isotopes and different final states in the hydrogen atom(s). In such calculations we found 
no limiting values for the momenta p e of photo-electron. In other words, all three-particle 
integrlas with spherical Bessel functions absolutely converge in this new method. Note also 
that in contrast with other similar methods in our approach we operate only with regular 
three-particle integrlas with Bessel functions, i.e. there is no need to discuss convergence of 
the Frullanian and/or truly singular integrals. 

This paper has the following structure. First, in the next Section we discuss the gen¬ 
eral formulas for the differential and total photodetachment cross-sections of the negatively 
charged hydrogen ion. It is shown that calculation of the photodetachment cross-section is 
reduced to numerical calculations of three-particle integrals with the spherical Bessel func¬ 
tions. Our approach applied in this study for calculations of such integrals is based on the 
Rayleigh’s formula for spherical Bessel function is discussed in Section III. Analytical for¬ 
mulas for the three-particle integrals with spherical Bessel functions are derived in Section 
IV. In Section V we discuss numerical computations of the photodetachment cross-sections 
of the negatively charged hydrogen ion(s). Concluding remarks can be found in the last 
Section. 


II. PHOTODETACHMENT CROSS-SECTION 

In [fj] we have derived the following formula for the differential cross-section of the pho¬ 
todetachment of the negatively charged hydrogen ion 


da 
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where /C is some numerical constant defined below, a — ^ — 7.2973525698 • 10~ 3 (~ -d^) 
is the dimensionless fine structure constant and «o ~ 0.52917721092 ■ 10 -8 cm is the Bohr 
radius (unless otherwise specified all numerical values of physical constants in this study are 


taken from 


12]). The k f vector is the unit vector which describes the photon propagation. 


Analogously, the unit vector n e = y- describe propagation of the free photo-electron which 
arises during photodetachment. Also, in this formula p e =| p e | is the momentum (absolute 
value) of the photo-electron and I n is the ionization potential of the negatively charged 
hydrogen ion determined in the case when the final hydrogen atom arises in the bound 
ns—state, where n > 1 and symbol s means that the angular momentum of this state equals 
zero, i.e. i = 0. The notation 7 Z^f in this formula is the radial part of the photodetachment 
amplitude. From this formula one finds the following formula for the photodetachment cross- 
section a ph of the negatively charged hydrogen ion(s): 
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For truly correlated three-particle wave functions the explicit formula for the radial part of 
the photodetachment amplitude is 

1 - N 
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3 0 30 J p 3 l —T32 | 

This amplitude is the overlap of the incident wave function of the H - ion (truly correlated) 
and the final state wave function which is the product of the radial part of the total wave 
function i? n o(r 3 i) of the final hydrogen atom H in one of its bound ns—states and the 
wave function of the moving photo-electron. For the ground state in the final hydrogen 
atom one finds 7?„o( r 3i) = 2exp(—r 3 i). If we can neglect by any interaction between this 
photo-electron and neutral hydrogen atom, then the wave function of the free photo-electron 
is a plane wave, which is represented as a combination of products of the spherical Bessel 

Is (see, e.g., Q, p.91). 


functions of the first kind with the corresponding Legendre polynomia 
This formula is called the Rayleigh expansion of the plane wave 


Conservation of the 
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angular momentum during photodetachment leads to the selection of only one £—component 
in such an expansion. 

It is shown below that with our choice of the incident and final wave functions the factor 
/C in Eqs.© - (ED equals ^. In this case the formulas Eqs.(J2D ~ © are written in the form 

1 — 20 ( Pe 
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These formulas are used in all calculations performed for this study. In the next two Sections 
we derive analytical formulas for the matrix element 7 Z^f which is included in the last two 
formulas. 


III. RAYLEIGH’S FORMULA FOR SPHERICAL BESSEL FUNCTIONS 


The method described below is, probably, the most effective and accurate approach which 
is used to determine the photodetachment cross-section of the negatively charged hydrogen 
ion(s). High efficiency and accuracy of this method directly follows from its simplicity and 
transparency. The method is based on the Rayleigh’s formula for the Bessel functions (uij ]. 
In the case of the ji(z) = jiipr^) Bessel function the Rayleigh’s formula takes the form 

djo(z 


3 i\z) = 


dz 


d /sin z 
dz ^ z 


For z = pr %2 one finds from this equation 

1 d 


jiiprm) = 


(jo(pr 


r 32 = -- 


1 r d / sin ( 7 / 1 - 32 ) 


(7) 


( 8 ) 


p 2 dr 32 v Ip 2 L<9r 32 v r 32 

where p = p e is the absolute value of momentum of the outgoing photo-electron p e , while r 32 
is the scalar distance between the hydrogen nucleus (heavy particle 3) and photo-electron 
(particle 2). The particle 1 designates another (first) electron which remains bound after 
photodetachment. The same system of notations is used everywhere below in this study. 

Let us designate the wave function of the incident H~ ion in the ground l 1 S —state 
by d'n-( r 32 , r 3 i> r 2 i), while the notation 0 n (^ 3 i) stands for the wave function of the final 
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ns —state of the hydrogen atom. Then we can write the following expression for the ampli¬ 
tude of the photodetachment jj5]: 


Mi^j = J (j) n {r 3 i)ji{pr 32 )^{r Z2 ,r 3ll r 21 )dV 


If r d -I 

= —y / ^n(r 3 i) yr— jo{pr 32 ) ^n-(r 32 ,r 31 ,r 2 i)dV 
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where dV = \r 2 1 dr-^dr^\dr 2 1 is an elementary volume in the relative coordinates r^ 2 , r 3 i 

and r 2 \. Each of the relative coordinates is the difference between the two Cartesian coor¬ 
dinates of the corresponding particles, e.g., r %3 =| iy — x 3 |. It follows from such a dehnition 
that: (a) relative coordinates are symmetric, i.e. r i3 = r 3 i , and (b) the following inequalities 
are always obeyed: | — r 3 k |< r l3 < rik + rjk for (■ i,j , k) = (1, 2, 3). Now, we can transform 

the expression in the right-hand side of Eq.(j9j) to the following form 


If 
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which contains the partial derivative of the incident wave function, i.e. the bound state wave 
function of the H~ ion. The function $( 732 , r 3 i, r 2 \) in Eq. flTO is related with the actual 
wave functions of the negatively charged hydrogen ion ^-(r^, r 3 i, ^ 21 ) by the following 
relation 


d 

®{r 32 ,r 3U r 21 ) = - —T H -(r 32 , r 31 , r 21 ) 
Or 32 


( 11 ) 


In this study the trial variational wave functions of the ground (singlet) l 1 *?—state of the 
negatively charged hydrogen ions are constructed in the following (exponential) form (see 
also 


1 ~ N 

T h - = -(1 + P 12 ) '52C i enp(-a i r 32 - ftr 3 i - ^r 21 ) 

z i =1 

1 , N 

= -(1 + P12) + A)%- {oti + li)u 2 - (^ + 7*)%)] (12) 

1 i =1 


where the notation P 12 stands for the permutation operator of identical particles (electrons), 
Ci {i = 1 , 2 ,..., N) are the linear parameters of the exponential expansion, Eq. (fT 2 lh while 
and 7 * are the non-linear parameters of this expansion. The trial wave function, 
Eq. (fT 2 lh has the correct permutation symmetry which is required for the spatial part of 
the singlet wave functions in the two-electron atomic systems. The non-linear parameters 
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in Eg. (112 j) must be varied in calculations to increase the overall efficiency and accuracy 
of the method. There are 3N additional conditions for these non-linear parameters are: 
cq + j3i > 0, oii + 7 j > 0, + 7 i > 0 (i = 1, 2,..., N). These conditions must always be 

obeyed to guarantee convergence of all three-particle integrals arising in computations. In 
this study we assume that all non-linear parameters in Eq. (1121) are the real numbers with 
zero complex parts. 

Note that Eq. (fl2l) is, in fact, the Laplace transform of the wave function T. It follows 
form Eq. (fl2]) that 
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The formulas, Eqs. (TT3|) - (TT4|) . are used in this study to determine the first and second partial 
derivatives of the trial wave functions. 


IV. ANALYTICAL FORMULAS FOR THREE-PARTICLE INTEGRALS WITH 
BESSEL FUNCTIONS 


In this Section we discuss analytical approach which was found to be very effective for 
analytical and numerical computations of the three-particle integral(s) derived above (see, 
Eq.((4])) which include one spherical Bessel function jo{pr 32 ). The three-particle integral with 
the spherical Bessel function jo(pr 32 ) which is needed in such calculations takes the form 




pZ 

1 


p 


exp(—c^r^ - A'ffil - 'yir 21 )j 0 (pr 32 )r 32 r 31 r 21 dr 32 dr 31 dr 21 


exp(—ctir 3 2 - Ann - 7 * r 2 i) sin (pr 32 )r 3 1 r 2 1 dr 32 dr 3 1 dr 2 1 (15) 


It is clear that analytical and/or numerical computations of three-particle integrals written in 
the relative coordinates is a difficult problem, since the three relative coordinates r 32 , r 31 and 
r 2 i are not truly independent. Indeed, three inequalities | Tik—Vjk |< r t] < r^+rjk are always 
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obeyed for three relative coordinates r 3 2 ,r 3 i and r 2 i (they are mentioned in the previous 
Section). In real three-body calculations it is better to use three perimetric coordinates 
u 3 ,u 2 ,u 3 which are simply related with the perimetric coordinates: w, = + r ik — 7 *+, 

where ( i,j,k ) = (1,2,3) and rjj = Vji. The inverse relations are: rij = Ui + Uj. Each of the 
three perimetric coordinates u t (i = 1,2,3) is: (1) non-negative, (2) truly independent from 
other perimetric coordinates, and (3) changes from zero to infinity. In perimetric coordinates 
the last integral from Eq. dlSp is written in the form 

2 f+oo r- 1-00 r+oo 

I(ai,(3i,”/i]p) = — / / exp [-(a* + (3i)u 3 - (a* + 7 +i 2 - (A + 7 i+i] x 

p z Jo Jo Jo 

sin (pu 2 + pu 3 )(u\ + uiu 2 + uiu 3 + u 2 u 3 )duidu 2 du 3 (16) 


where the factor 2 is the Jacobian of the transition from the relative to perimetric coordinates 
(t 32 , r 31 , r 2 i) —> (ui, u 2 , u 3 ). To derive the final expression we apply the known trigonometric 
formula sin (pu 2 + pu 3 ) = sin(pn 2 ) cos (pu 3 ) + cos (pu 2 ) sin(pn 3 ) and introduce the following 
notations: Z = oti + /?*, E = a* + 7 *, X — fa + 7 j. In these notations the integral Eq. (TT 6 ]) is 
written as the sum of eight integrals. 

To illustrate the process of integration in perimetric coordinates let us consider analytical 
computations of one of these eight integrals, e.g., 
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p 2 x 2 [y 2 + p 2 ] 2 z 2 + p 2 pX 2 [Y 2 +p 2 ] 2 (z 2 + p 2 ) 

Calculation of seven remaining integrals is also simple and transparent. Note that our 
analytical computations are based on the formulas from jl5|. The final expression for the 
integral /(«*, 7 *;p) takes the form 
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Note that the formula, Eq. (fl 8 |i . does not contain any singular and/or quasi-singular term 
represented by the Frullanian integrals jlfj, [tJ:] which must be regularized before actual 
numerical computations. This is an obvious advantage of our approach which is based on 
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the Rayleigh’s formula for spherical Bessel functions. It should be mentioned here that 
integrals of the general type / 0 °° R ax )~/( bx ) ( j x we re described by Italian mathematician G. 


Frullani in 1828 16]. These integrals look like singular, but they converge, if the finite values 
/(0)(^ 0) and /(oo) do exist and a and b are both positive (or negative) real numbers. In 
the photodetachment of the H~ ion one always finds three-particle Frullanian integrals which 
look like singular integrals, but can be reduced to a few infinite (but convergent!) sums of 
the regular three-particle integrals. 

In earlier studies H, H we have developed a number of different methods (six meth¬ 
ods) for analytical and/or numerical computations of three-particle integrals which contain 
Bessel functions of the first and second kind. Some of these methods are more successful 
in applications than others. In particular, we have found that methods based on power- 
series expansions of the spherical Bessel functions are not appropriate for highly accurate 
computations of the overlap integrals with the bound state wave functions. The reason is 
obvious, since contribution of the three-particle integrals rapidly increases when powers of 
the relative coordinate r %2 increase in such integrals. In general, any power-series expan¬ 
sion of the spherical Bessel functions contains large and very large powers of this relative 
coordinate r^ 2 - Convergence rate of such three-particle integrals slows down rapidly. For 
large n all integrals which contain r% 2 become divergent, if the contribution of high power 
of t 32 cannot be compensated by the increasing powers of some small parameter (s). This 
explains a restricted use of the power-series expansion for the spherical Bessel functions in 
actual computations. Our method developed in this study is free form these difficulties. 
Moreover, the formulas presented above are remarcably stable in calculations with large 
number of the non-linear parameters. This allows us to solve the problem to very high ac¬ 
curacy. Furthermore, by using these formulas we can also determine the photodetachment 
cross-section in those cases when the final hydrogen atom is formed in the excited states. 
The power-type series for Bessel functions derived in [18] and Q are used in calculations 
of the three-particle integrals with Bessel functions only to check our current formulas for 
small values of p e (< 0.015 — 0.030). 
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V. NUMERICAL CALCULATIONS 


Numerical calculations of the photodetachment cross-section of the negatively charged 
hydrogen ions H~ are performed in this study with the use of different variational wave 
functions, or trial wave functions, for short. These trial wave functions contain different 
number of the basis wave functions, Eq. (fT2]h i.e. they have different overall accuracy. In 
this study we have used trial wave functions, Eq. (fT2]h with N = 50, 100, 200, 300 and 350 
basis functions (exponents). Convergence of onr photodetachment cross-section upon the 
number of basis functions N in the trial wave functions is shown in Table I for the H~ ion 
with infinitely heavy nuclear mass which is designated as the °°H~ ion. 

Table II contains photodetachment cross-sections of the °°hU ion computed for 200 dif¬ 
ferent values of p = p e , which is the absolute value of momentum p e of the emitted photo- 

2 

electron. In our Tables the kinetic energy of the photo-electron E e = ^ is used, since the 
dependence a p (E e ) is directly measured in modern experiments. In Table II the kinetic 
energy of the photo-electron E e is related with the ‘index’ i as follows: 

r>? 

E e (i) = y , where p t = 0.005 + (i - 1) ■ 0.005 (19) 


in atomic units. The photodetachment cross-section of the H ion is written in the form: 
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where 77j_>/ is the radial part of the photodetachment amplitude, Eq.Q and I\ = E(°°H~) 
+ 0.5 a.u ., where E(°°hU) is the total energy of the H _ ion. The numerical value of the 
factor K in this formula equals ^2 which equals to the square of the norm of the angular 
part of the truly correlated, two-electron wave function of the m * ty H _ ion. Finally, one finds 
the following expression for the photodetachment cross-section 
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where the expression for 77is given by Eq.flTJ). The formula, Eq. fl2l]l . has been used in 
our calculations of the photodetachment cross-section of the H - ion (see Tables I and II). 
In general, our results from Tables I and II agree very well with the results of the earlier 
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and recent experimental studies 20 ] - 24], Agreement between our photodetachment cross- 
section of the H~ ion and results obtained in other theoretical/computational works (see, e.g., 


25]) can be considered as very good for those cases when the final hydrogen atom is formed 


in its ground Is—state. Note that experimental studies of the H ion photodetachment with 
the formation of the final hydrogen atom in the excited S'—states are extremely difficult to 


perform. Agreement of our photodetachment cross-sections determinec 


hydrogenic state and results known from earlier theoretical studies 


26 


in the case of 2 s 
and 27] is good. 


However, it is clear that to obtain a perfect agreement between theoretical cross-sections in 
this case one needs to conduct additional research and evaluations (see also discussion in 
the Conclusion). 

Another interesting question is to determine photodetachment cross-section of the H ion 
in the case when the final hydrogen atom arises in the first excited s— state. All calculations 
are almost identical to the case considered above when the final hydrogen atom is formed 
in its ground l 2 ^—state. However, in this case the ionization potential of the H ion is 
/ = i£(H - ) + 0.25 a.u. and in the formula for the photodetachment amplitude one needs to 
replace the bound state wave functions of the hydrogen atom, e.g., for the °°H atom: 


2 . exp(-r) -Lexp(-^)(l - Q (22) 

The photodetachment cross-section of the °°H~ ion in the case when the final hydrogen atom 
arises in the first excited s —state is determined with the same trial wave functions of the 
H~ ion (see above, for N = 350). The photodetachment cross-section (in cm 2 ) for this case 
is presented in Table III. 

Note that calculations of the photodetachment amplitude in the case of the first excited 
state (in the final hydrogen atom) requires computation of the two three-particle integrals 
with Bessel functions. One of these integrals essentially coincides with Eq. (fT5lh while another 
integral contains an extra power of r 3 i, i.e. 


J(a,/3,7;p) 


pZ 


exp(—cqr 32 - far 3i - 7i7"2i) sm(pr 3 2 )rl 1 r 2 idr 32 dr 31 dr 2 i (23) 


Analytical expression for this integral can be derived in the same way as we did to determine 
the integral /(a,/3, 7 ; p), Eq. (fT51h However, the final expression for the J(a,/ 3 , 7 ;p) integral 
is difficult, since it contains a large number of different terms. To determine the integral 
J(a,/3, 7 ;p) in this study we apply another method which is based on the use of numerical 
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differentiation of the integral I (a, (3,^]p), Eq. (fT5]) . upon the non-linear parameter a. To 
calculate this first order derivative we use the following formula 


d 

J(a,/3,7;p) = - — I {on, A, 7b p) 


12A L 


I {on - 2A, /3j, 7i;p) - 8 1{ai - 


+ 81(an + A, A,7bh0 - /(a* + 2A, A,7bP) 


(24) 


where the three-particle integral /(a,, fy, is defined by Eq. fflfUl . In Eq. (124)1 for each 
a.i the parameter A is defined as follows: A = 5 ■ ai, where i = 1,... ,N and 6 is a small 
number, e.g., 5 = 1 • ICE 6 , or 5 = 1 • ICE'. The overall accuracy of the formula, Eq.( 


can be evaluated as / • (see, e.g., 28]), where / is a positive numerical factor. To prove 


the correctness of numerical calculations of the first-order derivative one needs to perform 
a series of calculations by varying small parameter 6 from, e.g., 1 • ICE 4 up to 1 • ICE 8 (if 
calculations are performed with with the use of the quadruple precision). Such a series of 
calculations allows one to evaluate the total number of stable decimal digits in this derivative. 
In general, the computed derivative must be stable in 10 - 12 decimal digits (at least). 

As follows from comparison of Tables II and III the photodetachment cross-section for the 
ground state in the final hydrogen atom is significantly larger than analogous cross-section 
for the first excited state. This can be explained by ‘almost exact’ orthogonality of the 
two-electron wave function of the Hr ion \h H -( r 32 , Ei, r 2 i) and the 2s—hydrogenic orbital 
0(r 31 ). This fact is known for the natural orbital expansions constructed for the negatively 
charged atomic ions, e.g., for the hydrogen ion(s) H~. It is closely related to the fact that 
the negatively charged hydrogen ion has only one bound l 1 S— state. However, there is no 
such an ‘almost exact’ orthogonality between the T H -( r 32 ) r 3 i, r 2 i) function and 3s— and 
4s—hydrogenic orbitals. This means that the corresponding photodetachment cross-sections 
for the 3s— and 4s—states in the final hydrogen atom can be larger than for the 2s—state. 
Another interesting detail which has been noticed in our calculations is the presence of 
multiple maxima in photodetachment cross-sections of all excited ns —states for n > 3 (in 
the final hydrogen atom). Already, for the 3s—state in the hnal hydrogen atom one finds 
two such maxima in the photodetachment cross-section. This follows from the explicit form 
of the photodetachment amplitude 7Eq. (f4j) and the fact that the hydrogenic i?„o(r) 
functions (as well as spherical Bessel function jiijpef)) oscillate, i.e. vary between their 
maximal and minimal values for r > 0 (and n > 1). This interesting problem will be 
discussed in our next study. 
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The energy conservation law leads to the conclusion that formation of the final hydrogen 
atom in the excited atomic states, including excited s— states, requires substantially larger 
energy of the incident photon. However, if we assume that the actual stellar photosphere is at 
thermal equilibrium with T max fa 4000 - 8000 K , then the part of such ‘high-energy photons’ 
with energies > 80,000 K is extremely small. For instance, to produce photodetachment 
of the H - ion with the formation of the final hydrogen atom in the excited 2s—state the 
‘photon temperature’ Tf must be > 78,944 K. Analogous ‘photon temperature’ to form 
the final H atom in the excited 3s—state is around 118,416 K. For our Sun the fraction 
of such ‘hot’ thermal photons with T > 78944 K is very small (< 1 • ICC 2 %). Even in 
photospheres of late A—stars the fraction of ‘hot’ thermal photons is less than 0.5 %. In 
reality, such photons become noticeable in photospheres of B-stars with surface temperature 
T > 20, 000/\. However, in classical B-stars the photodetachment of the II ions does not 
play any role and most of the emitted radiation is absorbed by neutral hydrogen helium 
atoms and helium one-electron ions. This situation can be different in the Be-stars, i.e. 


EL 


in the B-stars which also include a rapidly rotating, low-dense disk of hydrogen atoms 
30]. Spatial radius of such a ‘hydrogen’ disk Rp is substantially larger than the radius of 
the central hot star R$. For a typical Be-star Rp fa 150 • 10 8 km , while R$ < 3 • 10 6 km. 
Illumination of this disk by radiation emitted by the central hot B-star produces excitations 
of hydrogen atoms which are moving in the stellar disk. Excitation of hydrogen atoms and 
following emission of absorbed radiation generates a very special optical spectra of the Be- 
stars, which contain thermal (continous) spectra of radiation emitted from the central hot 
B-star and a linear spectra which corresponds to the transition between excited states in the 
hydrogen atoms moving in the stellar disk. A well known Be-star is the 7 Cassiopeia which 
can easily be located and observed at the clear night sky in Northern Hemisphere. The 
emission spectrum of the 7 Cassiopeia contains a typical optical spectrum of a hot B-star 
and the red line which corresponds to the Balmer optimcal series in hydrogen atoms. Such 
a very intense red line essentially determines the spectrum of this star. In general, fa 35 % 
of all B-stars are the Be-stars, i.e. they have an additional low-dense disk which contains 
hydrogen atoms. For such stars photodetachment of the H~ ions with the formation of the 
final hydrogen atoms in their excited states may be important in some outer areas of stellar 
disk. In other stars the photodetachment of the H~ ions which leads to the formation of the 
final hydrogen atoms in the excited states is of only a restricted academic interest. 


13 



Finally, let us consider the photodetachment of the actual hydrogen ions with finite 
nuclear masses, i.e. ions of the protium 1 H~, deuterium 2 II and tritium 3 H . To consider 
photodetachment of these negatively charged ions one needs to take into account the effect 
of ‘nuclear recoil’, i.e. recoil of the central hydrogen atom. In this study we apply the same 
nuclear masses of the hydrogen isotopes which were used in j^]: 


M p = 


938.272046 

0.510998910 


1875.612859 _ 2808.920906 

0.510998910 ’ * “ 0.510998910 


(25) 


These nuclear masses are the most recent values obtained in high-energy experiments. The 
best-to-date total energies obtained for the negatively charged ions of the protium 1 H“, 
deuterium 2 H~ and tritium 41 with these nuclear masses can be found in Table I of jsj]. 
Numerical computations of the photodetachment cross-sections of the negatively charged 
hydrogen ions are more complicated than for the model °°H~ ion. General theory of atomic 
photodetachment leads to a conclusion that photodetachment cross-section of a less bounded 
electron is always smaller than analogous cross-sections of strongly bounded atomic elec¬ 
trons. On the other hand, the negatively charged hydrogen ions have the Is 2 —electron 
configuration. The less bounded electrons in this case mean larger ‘geometric radius’ of the 
negatively charged hydrogen ion, i.e. a larger cross-section of photodetachment. In reality, 
there is a competition of these two opposite factors in the following series of ions: °°H~, 3 H~ 
(tritium), 2 H~ (deuterium) and 1 H _ (protium). 

Let us discuss possible corrections to the photodetachment cross-section which are di¬ 
rectly related to the finite nuclear mass. There are two groups of such corrections. The 
first group contains the common factor 1 + jj (in atomic units), where M is the mass of 
the central nucleus. These corrections arise due to additional factor(s) in the atomic wave 
functions of the final hydrogen atoms. The second group of corrections includes the common 
factor 1 + jfyj- These corrections are directly related with the hydrogen atom recoil dur¬ 
ing photodetachment of the H _ ion. For the negatively charged hydrogen ions the energy 


conservation law takes the form: hu = 


_ Pe 


+ 


P 2 
r H 


+ /, where Ph =| P h | is the absolute 


2 m e 1 2(M+1) 

value of momentum of the neutral hydrogen atom and I is the ionization potential of the 
H~ ion. The conservation of the electron and nuclear momenta in the H~ ion leads to the 
relation: p e = —Ph, or p e = Ph- From these two equations one can write: 


hu = 


2 

Pe 


1 + 


2 m e L 2(M + 1)J 


+ / = 


■pi 
2 L 


1 + 


2(M + 1) 


+ 


7 


(26) 
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where 7 = y/2I and equality in the right-hand side of Eq. fl26l) is written in atomic units. 
This lead to the following expression for the differential cross-section of the photodetachment 
of the H~ ion [51: 


da,, = 


-aa n 


7r 


Pe 


Pet 1 + 


2(M+1) 




<hi(rai)ji(p e r32) 


] + 7 2 

,d_ _d_ 

^<9ri dr 


M + l 


x 


(27) 


) ^H-(r32,r3i,r2i)r32r3ir2idr32dr 3 idr2i(iO}| do 


where the factor (l + —hyj reflects the fact that the electron velocity in the hydrogen 

atom/ion with the finite nuclear mass is larger than in analogous hydrogen atom with 

infinitely heavy nucleus. There is also an additional mass dependent factor (l + which is 

included in the normalization constant of the bound state wave function of the final hydrogen 

atom <5 h(' 731 ) • In general, combination of a number of mass-dependent factors leads to 

relatively complex mass-dependence of the differential cross-section of the photodetachment 

of the ion. Results o four preliminary calculations of the photodetachment cross-sections 

of the 3 H~ (protium), 2 H _ (or D _ ) and 3 H~ (or T~) can be found in Table IV for different 

2 

energies of the photo-electron E e = ^ (in atomic units). In calculations of each ion shown 
in Table IV we have used 350 basis functions in Eq. (ll2|) . Photodetachment cross-section 
of the model ion °°H _ (determined with 350 basis functions) is also presented in Table IV. 
The corresponding total energies obtained with these trial wave functions are also listed in 
Table IV. Note that these energies are very close to the ‘exact’ vales obtained in [sj]. This 
indicates very high quality of the trial wave functions. As follows from Table IV these cross- 
sections are relatively close to the photodetachment cross-sections of the H~ ion (see Table 
II). However, some differences in such cross-sections do exist and must be investigated. We 
are planning to investigate such differences in our next study. 


VI. CONCLUSIONS 

We have considered the photodetachment of the negatively charged hydrogen ion(s). 
By using our new method we determine the photodetachment cross-sections of the H - 
ion with the infinitely heavy nucleus and for analogous ions of some hydrogen isotopes, 
including protium 1 H - , deuterium 2 H - (or D~) and tritium 3 H~ (or T _ ). The crucial 
part of such calculations is highly accurate computations of the three-particle integrals 
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which also contain spherical Bessel functions of the first kind and truly correlated wave 
functions of the two-electron hydrogen ions H~. In our method to solve all problems arising 
in this part we apply the Rayleigh’s formula for the spherical Bessel function j\ (peT'w)- 
Photodetachment cross-sections of the negatively charged °°H _ ions determined in this study 
perfectly coincide with the photodetachment cross-sections determined in experiments 20] 


- 24], Furthermore, the agreement between data from Table I in 24] and our cross-sections 
determined at the same energies (see Tables I and II) is close to absolute (or 1:1). In 
general, our photodetachment cross-sections computed in the ‘velocity representation’ are 


located between five experimental curves ploted in Fig.5 from 


24], In fact, our data are very 


close to the experimental curve ploted by small circles ‘o’, which correspond experimental 


data measured in 


24], It should be mentioned that the paper 24] was not known to me 


when I conducted these calculations (it was published after this paper has been submitted). 
Such an accurate coincidence indicates clearly the ‘velocity picture’ used by us and based 
on the fundamental QED rules is appropriate and accurate to consider photodetachment 
of the negatively charged hydrogen ion(s). Moreover, a very accurate coincidence of our 
(theoretical) and experimental data shows that overall contribution of the electron-electron 
correlations and other similar effects into the H~ photodetachemnt cross-sections is relatively 
small andean be neglected in the first approximation (as expected) (see discussions in JsJ and 


25]). In part, such an accurate agreement with the experimental data can be expained by the 


fact that in our calculations we have applied very accurate bound state wave functions of the 
H~ ion and its isotopes. For instance, the total energies of the °°H _ ion obtained with 300 
and 350 exponential basis functions are -0.527751016544180 a.u. and -0.527751016544278 
a.u. which coincide with the ‘exact’ energy (see, e.g., Q) in 13 decimal digits. 

In the case when the final hydrogen atom is formed in the excited 2 2 S —state we still 
do not have reliable experimental data, since any direct experimental study of the H~ ion 
photodetachment with the formation of the final hydrogen atom in the excited S'—states is 
extremely difficult. A few theoretical evaluations of the corresponding cross-section(s) { 25 1 
- [ 27 ]] were only approximate. In this study we consider photodetachment of the negatively 
charged hydrogen ions in those cases when the final hydrogen atom is formed in the excited 
2s—state. The corresponding cross-sections are determined with the use of numerical dif¬ 
ferentiation of the three-particle integral(s) in respect with one of its four parameters. The 
computed photodetachment cross-sections are in agreement with the cross-sections evalu- 
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ated earlier in 


26j and 


0 . 


However, we have to note again that the procedure of numerical 


differention was found to be sufhcintly accurate for our current purposes, but in the future 
highly accurate calculations it must be replaced by more accurate methods. Briefly, we can 
say that our current results from Table III can be considered as preliminary. 

Another problem discussed in this study is to determine photodetachment cross-sections 
of the negatively charged, two-electron ions of the hydrogen isotopes: 1 H _ (protium), 2 H~ 
(deuterium D~) and 3 H~ (tritium T _ ). In each of these ions the nuclear masses are finite. 
Recently, we have made a significant progress in study of the photodetachment of such ‘real’ 
hydrogen ions, but still a few additional steps are needed to obtain an absolute agreement 
with the known experimental data. In particular, we need to take into account all effect 
related with electron-electron correlations in the final state of the photodetachment process, 
Eq. (JT|) . In general, such corrections are small, but their overall contribution can noticeably 
change cross-sections at some energies. Our approach to include such corrections is described 


in detail in 


Analysis of such corrections will be considered in the next study. Another 


approach to this problem was discussed in [31 ]. 
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TABLE I: Photodetachment cross-section (in cm 2 ) of the negatively charged °°H~ ion in the ground 
l 1 S— state. N is the total number of basis function in the trial wave function and E is the total 
energy (in a.u .) for this wave function. The final hydrogen atom is formed in the ground ls-state. 
The notation E e stands for the energy of the emitted photo-electron (in a.u.). 


Eh- 

-0.5277751015678648 

-0.5277751016514754 

E e 

*{E e )(N = 50) 

a(E e )(N = 100) 

0.0000125 

0.29605783690E-17 

0.29598715990E-17 

0.0001125 

0.87973761028E-17 

0.87953808289E-17 

0.0003125 

0.14386547859E-16 

0.14383603423E-16 

0.0006125 

0.19580556054E-16 

0.19577122883E-16 

0.0010125 

0.24256110179E-16 

0.24252654008E-16 

0.0015125 

0.28320090502E-16 

0.28316982705E-16 

0.0021125 

0.31712467227E-16 

0.31709927491E-16 

0.0028125 

0.344062.50159E-16 

0.34404338215E-16 

0.0036125 

0.36404882768E-16 

0.36403529204E-16 

0.0045125 

0.37737746241E-16 

0.37736804077E-16 

0.0055125 

0.38454542597E-16 

0.38453841338E-16 

0.0066125 

0.38619295849E-16 

0.38618684920E-16 

0.0078125 

0.38304584444E-16 

0.38303959740E-16 

0.0091125 

0.37586442834E-16 

0.37585756010E-16 

0.0105125 

0.36540186616E-16 

0.36539440563E-16 

0.0120125 

0.35237254090E-16 

0.35236489739E-16 

0.0136125 

0.33743033576E-16 

0.33742313147E-16 

0.0153125 

0.32115564584E-16 

0.32114955391E-16 

0.0171125 

0.30404958234E-16 

0.30404.519702E-16 

0.0190125 

0.28653369842E-16 

0.28653144950E-16 

0.0210125 

0.26895364963E-16 

0.26895376194E-16 

0.0231125 

0.25158540750E-16 

0.25158789884E-16 

0.0253125 

0.23464290471E-16 

0.23464761154E-16 

0.0276125 

0.21828625603E-16 

0.21829287584E-16 

0.0300125 

0.2026299421.5E-16 

0.20263808223E-16 

0.0325125 

0.18775054831E-16 

0.18775977338E-16 

0.0351125 

0.17369381282E-16 

0.17370368625E-16 

0.0378125 

0.16048086277E-16 

0.16049097863E-16 

0.0406125 

0.14811360133E-16 

0.14812360646E-16 

0.0435125 

0.13657926934E-16 

0.13658887600E-16 

0.0465125 

0.12.585423999E-16 

0.12586323045E-16 

0.0496125 

0.11590712460E-16 

0.11591534944E-16 

0.0528125 

0.10670127572E-16 

0.10670864770E-16 

0.0561125 

0.98196773647E-17 

0.98203258858E-17 

0.0595125 

0.90351977918E-17 

0.90357585655E-17 

0.0630125 

0.83124717978E-17 

0.83129490365E-17 

0.0666125 

0.76473188426E-17 

0.76477190671E-17 

0.0703125 

0.70356605505E-17 

0.70359917250E-17 

0 07411 25 

0 64735672865E-17 

0 64738380909E-1 7 


-0.527751016542222 


o(E e )(N = 200) 
0.29-596998756E-17 
0.87949302746E-17 
0.14382993623E-16 
0.19576430332E-16 
0.24251915406E-16 
0.28316242717E-16 
0.31709256092E-16 
0.34403814787E-16 
0.36403217986E-16 
0.37736736862E-16 
0.38454011788E-16 
0.38619053036E-16 
0.3830446282.5E-16 
0.37586322075E-16 
0.36540000770E-16 
0.35236987673E-16 
0.33742710013E-16 
0.32115231333E-16 
0.30404672073E-16 
0.28653184692E-16 
0.26895323402E-16 
0.25158669419E-16 
0.23464598890E-16 
0.21829107569E-16 
0.20263630834E-16 
0.18775818331E-16 
0.17370238921E-16 
0.16049003856E-16 
0.14812304846E-16 
0.13658869441E-16 
0.12.586339722E-16 
0.11591582203E-16 
0.10670937586E-16 
0.98204190031E-17 
0.90358669036E-17 
0.83130679640E-17 
20G478445644E-17 
0.70361204564E-17 
0 64739674087E-1 7 


-0.527751016544180 


a(E e )(N = 300) 
0.29598550732E-17 
0.87953448276E-17 
0.14383551527E-16 
0.19577016434E-16 
0.24252450461E-16 
0.28316688061E-16 
0.31709607249E-16 
0.34404089266E-16 
0.36403439874E-16 
0.37736924254E-16 
0.38454171384E-16 
0.38619181632E-16 
0.38304552.518E-16 
0.37586365790E-16 
0.36539996014E-16 
0.35236937881E-16 
0.33742623905E-16 
0.32115121065E-16 
0.30404551036E-16 
0.28653065525E-16 
0.26895216560E-16 
0.25158582386E-16 
0.23464535980E-16 
0.21829070170E-16 
0.20263617929E-16 
0.18775827151E-16 
0.17370265606E-16 
0.16049044055E-16 
0.14812354194E-16 
0.13658923911E-16 
0.12586395851E-16 
0.11591637205E-16 
0.10670989381E-16 
0.98204661807E-17 
0.90359086434E-17 
0.83131039366E-17 
0.76478748231E-17 
0.70361453329E-17 
0 61739874186F,-17 


-0.527751016544278 


a(E e )(N = 350) 
0.29598566799E-17 
0.87953664283E-17 
0.14383597731E-16 
0.19577077184E-16 
0.242.52520229E-16 
0.28316758397E-16 
0.31709665345E-16 
0.3440412.5710E-16 
0.36403454405E-16 
0.3773692.5369E-16 
0.384.54171955E-16 
0.38619193960E-16 
0.38304585002E-16 
0.37586421743E-16 
0.36540074168E-16 
0.35237033798E-16 
0.33742731565E-16 
0.32115234221E-16 
0.30404664124E-16 
0.28653174140E-16 
0.26895317595E-16 
0.25158673936E-16 
0.23464617137E-16 
0.21829140775E-16 
0.20263678346E-16 
0.18775878074E-16 
0.17370307915E-16 
0.16049078706E-16 
0.14812382158E-16 
0.13658946126E-16 
0.12586413199E-16 
0.11591650497E-16 
0.10670999349E-16 
0.98204734767E-17 
0.90359138388E-17 
0.83131075229E-17 
0.76478772150E-17 
0.70361468723E-17 
0 64739883797E-1 7 










TABLE II: Photodetachment cross-section (in cm 2 ) of the negatively charged hydrogen ion °°H 
in its ground l 1 ^—states. The final hydrogen atom is formed in the ground s—state. 


















TABLE III: Photodetachment cross-section (in cm 2 ) of the negatively charged hydrogen ion °°H 
in its ground l 1 ^—states. The final hydrogen atom is formed in the first excited -s—state. 


Pe 

a{E e ) 

Pe 

a{E e ) 

Pe 

a(E e ) 

Pe 

a(E e ) 

0.01 

0.15169453769E-18 

0.16 

0.16689095636E-17 

0.31 

0.21171774945E-17 

0.46 

0.18258475296E-17 

0.02 

0.29513081491E-18 

0.17 

0.17180044561E-17 

0.32 

0.21231650513E-17 

0.47 

0.17848013982E-17 

0.03 

0.42950207831E-18 

0.18 

0.17632152811E-17 

0.33 

0.21252205995E-17 

0.48 

0.17424897768E-17 

0.04 

0.55759395646E-18 

0.19 

0.18052224207E-17 

0.34 

0.21232787756E-17 

0.49 

0.16991672981E-17 

0.05 

0.68210106756E-18 

0.20 

0.18445390636E-17 

0.35 

0.21173257734E-17 

0.50 

0.16550750965E-17 

0.06 

0.80402375798E-18 

0.21 

0.18815119924E-17 

0.36 

0.21073985814E-17 

0.51 

0.16104391465E-17 

0.07 

0.92268403351E-18 

0.22 

0.19163339790E-17 

0.37 

0.20935824448E-17 

0.52 

0.15654691307E-17 

0.08 

0.10364910897E-17 

0.23 

0.19490634468E-17 

0.38 

0.20760069165E-17 

0.53 

0.15203577762E-17 

0.09 

0.11437296099E-17 

0.24 

0.19796480642E-17 

0.39 

0.20548408558E-17 

0.54 

0.14752805973E-17 

0.10 

0.12430738115E-17 

0.25 

0.20079498256E-17 

0.40 

0.20302867197E-17 

0.55 

0.14303959796E-17 

0.11 

0.13338055308E-17 

0.26 

0.20337698926E-17 

0.41 

0.20025744610E-17 

0.56 

0.13858455425E-17 

0.12 

0.14158199181E-17 

0.27 

0.20568720138E-17 

0.42 

0.19719553114E-17 

0.57 

0.13417547244E-17 

0.13 

0.14895139123E-17 

0.28 

0.20770037429E-17 

0.43 

0.19386956840E-17 

0.58 

0.12982335353E-17 

0.14 

0.15556298503E-17 

0.29 

0.20939149690E-17 

0.44 

0.19030713856E-17 

0.59 

0.12553774311E-17 

0.15 

0.16150991726E-17 

0.30 

0.21073734863E-17 

0.45 

0.18653622838E-17 

0.60 

0.12132682685E-17 
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TABLE IV: Photodetachment cross-section (in cm 2 ) of the negatively charged hydrogen ion in 
their ground l 1 ^—states. N is the total number of basis function in the trial wave function and 
E is the total energy (in a.u.) for this wave function. The final hydrogen atom is formed in the 
ground ls-state. The notation E e stands for the kinetic energy of the emitted photo-electron (in 
a.u.). 



°°H- 

3 h- 

2 H“ 

x h- 

Eh- 

-0.5277751016544278 

-0.527649048201840 

-0.527598324689626 

-0.527445881119668 

E e 

a{E e )(N = 350) 

cr(E e )(N = 350) 

cr(E e )(N = 350) 

a(E e )(N = 350) 

0.0000500 

0.2959856679870E-17 

0.5898071851203E-17 

0.5897932817961E-17 

0.58970881335E-17 

0.0001125 

0.5898578992686E-17 

0.8794584542552E-17 

0.8794335668102E-17 

0.87930359069E-17 

0.0002000 

0.8795366428292E-17 

0.1162901145802E-16 

0.1162861725426E-16 

0.11626832132E-16 

0.0003125 

0.1163008685588E-16 

0.1438220490456E-16 

0.1438163186360E-16 

0.14379323906E-16 

0.0004500 

0.1438359773146E-16 

0.1703629503788E-16 

0.1703551495594E-16 

0.17032638389E-16 

0.0006125 

0.1703803734141E-16 

0.1957494491642E-16 

0.1957393521732E-16 

0.19570436466E-16 

0.0008000 

0.1957707718383E-16 

0.2198356259143E-16 

0.2198230472500E-16 

0.21978124069E-16 

0.0010125 

0.2198613138832E-16 

0.2424946573180E-16 

0.2424794299533E-16 

0.24243016855E-16 

0.0012500 

0.2425252022884E-16 

0.2636199692909E-16 

0.2636019272279E-16 

0.26354456186E-16 

0.0015125 

0.2636558567114E-16 

0.2831258986430E-16 

0.2831048683222E-16 

0.28303876164E-16 

0.0018000 

0.2831675839673E-16 

0.3009478807359E-16 

0.3009236812510E-16 

0.30084823091E-16 

0.0021125 

0.3009957746377E-16 

0.3170421932053E-16 

0.3170146428578E-16 

0.31692930205E-16 

0.0024500 

0.3170966534484E-16 

0.3313852966036E-16 

0.3313542227316E-16 

0.33125851743E-16 

0.0028125 

0.3314466238433E-16 

0.3439728216636E-16 

0.3439380711190E-16 

0.34383161403E-16 

0.0032000 

0.3440412571010E-16 

0.3548182594435E-16 

0.3547797079138E-16 

0.35466220892E-16 

0.0036125 

0.3548939831482E-16 

0.3639514145850E-16 

0.3639089737677E-16 

0.36378024707E-16 

0.0040500 

0.3640345440464E-16 

0.3714166832346E-16 

0.3713703055060E-16 

0.37123027354E-16 

0.0045125 

0.3715072721740E-16 

0.3772712159963E-16 

0.3772208966385E-16 

0.37706959091E-16 

0.0050000 

0.3773692536921E-16 

0.3815830229366E-16 

0.3815288002287E-16 

0.38136635946E-16 

0.0055125 

0.3816884343565E-16 

0.3844290725784E-16 

0.3843710261617E-16 

0.38419769197E-16 

0.0060500 

0.3845417195483E-16 

0.3858934304742E-16 

0.3858316784240E-16 

0.38564778897E-16 

0.0066125 

0.3860131140057E-16 

0.3860654758202E-16 

0.3860001708227E-16 

0.38580615273E-16 

0.0072000 

0.3861919396038E-16 

0.3850382270889E-16 

0.3849695520654E-16 

0.38476591109E-16 

0.0078125 

0.3851711620594E-16 

0.3829068002096E-16 

0.3828349636687E-16 

0.38262227446E-16 

0.0084500 

0.3830458500200E-16 

0.3797670157192E-16 

0.3796922470517E-16 

0.37947114237E-16 

0.0091125 

0.3799117829127E-16 

0.3757141647769E-16 

0.3756367096565E-16 

0.37540786939E-16 

0.0098000 

0.3758642174276E-16 

0.3708419381512E-16 

0.3707620541543E-16 

0.37052619444E-16 

0.0105125 

0.3709968167445E-16 

0.3652415173423E-16 

0.3651594698662E-16 

0.36491733251E-16 

0.0112500 

0.3654007416812E-16 

0.3590008229333E-16 

0.3589168814598E-16 

0.35866922392E-16 

0.0120125 

0.3591638988830E-16 

0.3522039120618E-16 

0.3521183467941E-16 

0.35186593282E-16 
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